An Exact SU(2) Symmetry and Persistent Spin Helix in a Spin-Orbit Coupled System 
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Spin-orbit coupled systems generally break the spin rotation symmetry. However, for a model with 
equal Rashba and Dresselhauss coupling constant (the ReD model), and for the [110] Dresselhauss 
model, a new type of 517(2) spin rotation symmetry is discovered. This symmetry is robust against 
spin-independent disorder and interactions, and is generated by operators whose wavevector depends 
on the coupling strength. It renders the spin lifetime infinite at this wavevector, giving rise to a 
Persistent Spin Helix (PSH). We obtain the spin fluctuation dynamics at, and away, from the 
symmetry point, and suggest experiments to observe the PSH. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15. Gd 



The physics of systems with spin-orbit coupling has 
generated great interest from both academic and prac- 
tical perspectives 1 ". In particular, spin-orbit coupling 
allows for purely electric manipulation of the electron 
spii£iy&, an( j comc i De f use m ar eas from spintronics 
to quantum computing. Theoretically, spin-orbit cou- 
pling reveals fundamental physics related to topological 
phases and their applications to the intrinsic and quan- 
tum spin Hall effect 6 i 7 i 8 i 9 i 10 . 

While strong spin-orbit interaction is useful for ma- 
nipulating the electron spin by electric fields, it is also 
known to have the undesired effect of causing spin- 
decoherenc o 11 ! 12 . The decay of spin polarization reflects 
the nonconservation of the total spin operator, S, i.e. 
[S, 7i] 7^ 0, where Ti is any Hamiltonian that contains 
spin-orbit coupling. In this Letter, we identify an ex- 
act SU (2) symmetry in a class of spin-orbit systems that 
renders the spin lifetime infinite. The symmetry involves 
components of spin at a finite wave vector, and is different 
from the U(l) symmetry that has previously been asso- 
ciated with this class of modelsi 3 ^. As a result of this 
symmetry, spin polarization excited at a certain "magic" 
wavevector will persist. If this symmetry can be realized 
experimentally, it may be possible to manipulate spins 
through spin-orbit coupling without spin-polarization de- 
cay, at length scales characteristic of today's semiconduc- 
tor lithography. 

We consider a two-dimensional electron gas without 
inversion symmetry, allowing spin-orbit coupling that is 
linear in the electron wavevector. The most general form 
of linear coupling includes both Rashba and Dresselhaus 
contributions: 



H = h a(k y a x - k x a y ) 

2m 



kyCFy), (1) 



where k x ,y is the electron momentum along the [100] and 
[010] directions respectively, a, and j3 are the strengths 
of the Rashba, and Dresselhauss spin-orbit couplings and 
m is the effective electron mass. We shall be inter- 
ested in the special case of a = (3, which may be ex- 



perimentally accessible through tuning of the Rashba 
coupling via externally applied electric fields^. When 
a = /3, the spin-orbit coupling part of the Hamilto- 
nian reads a(k x + k y )(a x — <j y ). Rotating the spa- 
tial coordinates by introducing k± = -^(k y ± k x ), and 
also performing the global spin rotation generated by 
U = 77|(1 + -j7| (<t x + brings the Hamiltonian to 

the diagonal form: 



ft ReD = U^HU 



2m 



2ak 4 



(2) 



Throughout this paper, we shall be using both the orig- 
inal spin basis and the transformed spin basis for the 
7iR e D model, depending on the context. It should al- 
ways be remembered that o~ z in the transformed spin ba- 
sis corresponds to -^={<r x — a y ) in the original spin basis. 
^ReD is mathematically equivalent to the Dresselhauss 
[110] model, describing quantum wells grown along the 
[110] direction, whose Hamiltonian is given by: 



-\- k y 

n[u ] = — ^ 2ak x a z . 



(3) 



As the Hamiltonians 7i[R C D] and W[no] are already di- 
agonal, the energy bands (for 7i[R C D]) are simply given 
by: 



2m 



± 2ak +1 



(4) 



where [,] are the spin components in the new, unitary 
transformed, spin basis: S -> WSU. The bands in Eq.g] 
have an important shifting property: 



e l (k)=e r (k + Q), 



(5) 



where Q+ = 4ma, Q_ = for the 7Y[r c d] model and 
Q x = Ama, Q y = for the 7i[no] model. The Fermi sur- 
faces consist of two circles shifted by the "magic" shifting 
vector Q, as shown in Fig.pQ. 
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FIG. 1: [Color Online] Fermi surfaces of the model consists 
of two circles shifted by the wavevector Q — {Q+,Q~) = 
(4ma, 0). 

W[R e D] and 7i[no] have the previously known U(l) 
symmetry generated by a^2^, which ensures long life- 
time for uniform spin polarization along the z-axisi^. The 
exact SU (2) symmetry discovered in this work is gener- 
ated by the following operators, expressed in the trans- 
formed spin basis as: 

S Q = Efe C ll C k+QV S Q = °\+Q A C k{ 

S ° = ^k c lfk] ~ % c kl> ( 6 ) 

with Ck]\i being the annihilation operators of spin-up and 
down particles. These operators obey the commutation 
relations for angular momentum, 

[Sq,Sq] = ±.2Sq\ [Sq,Sq] — Sq (7) 

The shifting property Eq. ensures that the operators 
defined in Eq.j^] commute with the Hamiltonian, 

[WaaD, c^cgj = (e T (fc + Q)-q(fc))4 + Q T c £| =° ( 8 ) 

and similarly for ct cg^^, thus uncovering the SU(2) 

symmetry. This symmetry is robust against both spin- 
independent disorder and Coulomb (or other many-body) 
interactions as the spin operators commute with the 
finite- wave vector particle density p q — J^k c \ + - c k : 

[p q ,S±} = [p q ,SZ} = 0. (9) 

As a result single-particle potential scattering terms like 
^2 q V q p q , as well as many-body interaction terms like 
Tliq VqPqPqi a h commute with the generators ©, and the 
SU(2) symmetry is robust against these interactions. 

The SU(2) conservation laws imply that the expec- 
tation values of S z , Sq, and Sq have infinite lifetime. 
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FIG. 2: [Color Online] Spin configurations for two particles A 
and B depend only on the distance traveled along the x+ axis 
and not on their initial momenta. For the same x+ distance 
traveled, the spin precesses by exactly the same angle. After 
a length L+ = 2nh/Q+ the spins all return exactly to the 
original configuration. 



Since S z is a conserved quantity, a fluctuation of the z- 
component of spin polarization with wavevector q can 
only decay by diffusion, which takes time r s = 1/D s q 2 , 
where D s is the spin diffusion constant. Notice that the 
hermitian operators S X (Q) — ^(Sq + Sq) and S y (Q) = 

^(Sq — Sq) create a helical spin density wave in which 
the direction of the spin polarization rotates in the x, y 
plane in the transformed spin basis. Another manifesta- 
tion of this symmetry is the vanishing of spin-dependent 
quantum interference (or weak anti-localization) at the 
SU{2) poinOa. 

The infinite lifetime of the S x , S y spin helix with wave 
vector Q results from the combined effects of diffusion 
and precession in the spin-orbit effective field. A physical 
picture for the origin of the persistent spin helix (PSH) 
is sketched in Fig.j2], expressed in the transformed spin 
basis. Consider a particle propagating in the plane with 
momentum k. In time t, it travels a distance along the 
x + direction L + = k+t/m, while its spin precesses about 
the z axis by an angle 9 = Aak+t/h. Eliminating t from 
the latter equation yields 9 = 4amL + / H, demonstrating 
that the net spin precession in the x, y plane depends 
only on the net displacement in the x + direction and is 
independent of any other property of the electron's path. 
Electrons starting with parallel spin orientation and the 
same value of x+ will return exactly to the original ori- 
entation after propagating L + = 2nh/Q + . 

For sake of clarity, we have depicted the PSH in Fig. [3] 
in the original basis. For a range of values in different 
materials from weak to strong spin-orbit splitting a = 
10 2 — * 10 4 (m/s), the characteristic wavelength of the 
PSH is 10/im — > lOOnm. In GaAs, the typical value is 
L + = lfim. This characteristic wave length is on the 
scale of typical lithographic dimensions used in today's 
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semiconductor industry. 
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FIG. 3: [Color Online] (a) PSH for the HrcD model. The 
spin-orbit magnetic field is in-plane (blue), whereas the spin 
helix, for the choice of the relative signs in Eq.0 is in the 
(x+,z) plane, (b) PSH for the Timo] model. The spin-orbit 
magnetic field B^.' b " t , in blue, is out of plane, whereas the spin 
helix, in red, is in-plane. 

Mathematically, the PSH is a direct manifestation of a 
"non-abelian" flux in the ground state of the 7iR e D and 
the H[no] models. Following Ref. |r7)h we express WReD 
in the form of a background non-abelian gauge potential, 



WrteD = t; 1" t; — (M+ ~ 2maa z ) 2 + const. (10) 

2m 2m 

In contrast to the general case where the non-abelian 
gauge potential leads to a finite non-abelian field 
strength, the field strength vanishes identically for 
a = (3. Therefore, we can eliminate the vec- 
tor potential by a non-abelian gauge transformation: 
"ff £_) — > e^q>(i2max + )^^(x + , X-), $>i(x+,x—) — > 
exp(— i2max+) i &i(x+, Under this transformation, 

the spin-orbit coupled Hamiltonian is mapped to that of 
the free Fermi gas. The cost of the transformation is 
two-fold. First, the new wave function ty(x + ,x-) satis- 
fies twisted spin boundary conditions 18 . However, these 
have no physical effect, since a metallic system without 
off-diagonal- long-range-order is insensitive to a change of 
boundary conditions in the thermodynamic limit. Sec- 
ond, while diagonal operators such as the charge n 
and S z remain unchanged, off-diagonal operators, such 
as S~(x) — ip\(S)ipf(x) an d S + (x) — ip^(x)ipi(x) are 
transformed: S~(x) — > exp(— i4max+)S~ (x), S + (x) — > 
exp(iAmax+)S + (x). In the transformed basis, all three 
components of the spin obey the simple diffusion equa- 
tion, as the Hamiltonian is just that for free electrons 
without spin-orbit coupling. Hence S x = const is a 
solution to the diffusion equation. Transforming back 



demonstrates that this solution corresponds to (S x , S y ) — 
(cos(4maa;+), sin(4mair + )) = const ., which is the PSH. 

While spin is conserved at the SU{2) point, dynamics 
emerge when the conditions q = Q and/or a = (3 are 
not met. Solving for the coupled charge and spin dy- 
namics requires Boltzman transport equations, which we 
obtain below using the Kcyldish formalismiS^. Assum- 
ing isotropic scattering with momentum lifetime r, the 
retarded and advanced Green's functions are: 

G R > A (k,e) = (e-H±±)- 1 - (H) 

We introduce a momentum, energy, and position depen- 
dent charge-spin density which is a 2 x 2 matrix g(k, r, t). 
Summing over momentum: 

pM = / (0^ (fc ' r ' t} ' (12) 

gives the real-space spin-charge density p(r, t) = n(r, t) + 
S 1 (r,t)o-i, where n(r,i) and S*(r,t) are the charge and 
spin density and v = m/2ir is the density of states in 
two-dimensions. p{r,t) and g(k,r,t) satisfy a Boltzman- 
type equatio n 19 ; 20 : 

(13) 

that we now solve for the Hamiltonian of Eq.pQ for ar- 
bitrary a, p. To obtain the spin-charge transport equa- 
tions, we follow the sequence: time-Fourier transform the 
above equation; find a general solution for g(k,r,t) in- 
volving p(r, t) and the fc-dependent spin-orbit coupling; 
perform a gradient expansion of that solution (assuming 
d r « kp where kp is the Fermi wavevector) to sec- 
ond order; and, finally, integrate over the momentum. 
The anisotropic nature of the Fermi surfaces when both 
a and (3 are nonzero introduces coupling of charge and 
spin degrees of freedom beyond what was found with only 
Rashba coupling 21 . The final result, expressed in the 
i = (x+,X-,z) spatial coordinates and the original spin 
basis, is: 

8 t n = Ddfn + B 1 d x+ S x _ - B 2 d x _S x+ (14) 

d t S x _ = D8fS x _ + B x d x +n- C 2 d x _S z ~ T 2 S X _ (15) 

d t S x+ = DdfS x+ - B 2 d x _n - C x d x+ S z - T,S X+ (16) 

d t S z = DdjS z + C 2 d x _S x _ + dd x+ S x+ - (Ti + T 2 )S Z . 

(17) 

With an effective kp defined as \ / 2mEp, the constants 
in the above equations are: 

Bi = 2(a- f3) 2 (a + f3)kpr 2 , 

B 2 = 2{a + {3) 2 (a-f3)k 2 F T 2 , 

C 1 =2i(a + /3)fc|r, C 2 = 2±(a - (3)k F T 

T 1 = 2(a + (3) 2 k 2 F T, T 2 = 2(a-(3fk F T. (18) 



4 



The diffusion constant D = Vpr/2, Our results for 
the coupling coefficients are valid in the diffusive limit 
akpT << 1, (3kFi~ << 1 and reduce to the appropriate 
limits in the cases (3 — > or a — * O^iSi. We observe that 
for a general direction of propagation in the x± plane, the 
three components of spin, and the charge, are all coupled. 

Motivated by experiments that probe spin dynamics 
optically, we focus on the behavior of the out-of-plane 
component of the spin, which is S z in the original ba- 
sis. We assume a space-time dependence proportional to 
exp[i(ujt — q ■ r)], and compare q parallel to [110] and q 
parallel to [110]. 

The spin polarization lifetime for q parallel to [110] 
is not enhanced, as no shifting property exists in this 
direction. For this orientation, the four equations [14-17] 
separate into two coupled pairs, with n being coupled 
to S x , and S z coupled to S x _. The optical probe will 
detect the characteristic frequencies that are solutions to 
the equation that couples S z with S x _ : 



Transient spin-grating experiment, 



,23,24 



are particu- 



-Dq 2 



T^Jlf + WC*). (19) 



At a = (3 we have T 2 = C2 = and the characteristic 
frequencies become iu>i = —Dq 2 — T\, iu>2 = —Dq 2 . 

The spin polarization dynamics are qualitatively dif- 
ferent in the [110] direction, which is the direction along 
which the Fermi surfaces are shifted. In this case the 
four equations again decouple into two pairs, with n cou- 
pled to S x _ and S z coupled to S x+ . The characteristic 
frequencies are: 



1 



«Ji, a = -Dq' - |(2Tj + T 2 ± ^T* + 4q 2 Cl). (20) 

We note that iu>2 (q) has a minimum at a wavevector that 
depends only on a,f3, and m, i.e., is independent of fej? 
and r. When a = j3 the decay rates become iwi t 2 = 
— Dq 2 — Ti^fCiq, and iu>2(o) is zero (corresponding to the 
PSH) at the shifting wave-vector 4ma. We have checked 
that this relation continues to hold as we include higher 
order corrections to the transport coefficients. Earlier 
calculations based on the pure Rashba modelSiiSS also 
predict an increased spin lifetime at a finite wave vector 
(•v/l5/2)cmi, but the lifetime is enhanced relative to q = 
only by the factor 16/7. 



larly well-suited to testing our theoretical predictions 
(Eqs. [19,20]) and discovering the PSH, as they inject fi- 
nite wavevector spin distributions. If at t = an S z po- 
larization proportional to cos {qx+) is excited, we predict 
a time evolution, 



S z (q,t) = Ai{q)e iUl ^ t + A 2 {q)e i ^ t , (21) 



where 



^1,2(0) = \ 



1± 



(22) 



4q 2 C 2 

Thus, according to theory, the photoinjected spin po- 
larization wave will decay as a double exponential, with 
characteristic decay rates uj\^{q). As a — > (3 the weight 
factors of the two exponentials become equal and the de- 
cay rate of the slower component — * at the "magic" 
wave-vector. If transient grating measurements verify 
these predictions, they will be able to provide rapid and 
accurate determination of the spin-orbit Hamiltonian, en- 
abling tuning of sample parameters to achieve a — [3. 

In conclusion, we have discovered a new type of spin 
SU(2) symmetry in a class of spin-orbit coupled mod- 
els including the ReD model and the Dresselhauss [110] 
model. Based on this symmetry, we predict the exis- 
tence of a Persistent Spin Helix. The lifetime of the PSH 
is infinite within these models, but of course there will 
be other relaxation mechanisms which would lead to its 
eventual decay. We obtained the transport equations for 
arbitrary strength of Rashba and Dresselhauss couplings 
and these equations provide the basis to analyze experi- 
ments in search of the PSH. 
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